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ABSTRACT 

We have investigated the dynamic response of a coronal heimet streamer to the emergence from below 
of a current with its magnetic field in a direction opposite to the overlying streamer field. Once the 
emerging current moves into the closed region of the streamer, a current sheet forms between the emerg- 
ing field and the streamer field, because the preexisting field and the newly emerging field have opposite 
polarities. Thus magnetic reconnection will occur at the flanks of the emerged structure where the 
cu " c « l densit ? * s maximum. If the emerging current is large enough, the energy contained in the current 
and the reconnection will promptly disrupt the streamer. If the emerging current is small, the streamer 
will experience a stage of slow evolution. In this stage, slow magnetic reconnection occurring at the 
flanks of the emerged structure leads to the degeneration of the emerged current to a neutral point. 
Above this point, a new magnetic bubble will form. The resulting configuration resembles an mverse- 
po anty prominence. Depending on the initial input energy of the current, the resulting structure will 
either remain in situ, forming a quasi-static structure, or move upward, forming a coronal transient 
similar to coronal jets. The numerical method used in this paper can be used to construct helmet 
streamers containing a detached magnetic structure in their closed field region. The quasi-static solution 
may serve as a preevent corona for studying coronal mass ejection initiation. 

Subject headings : MHD — Sun: corona — Sun: prominences 


1. INTRODUCTION 

Considerable interest and work have recently been 
devoted to studying the complex interaction of coronal 
helmet streamers, coronal mass ejections (CMEs), and 
erupting prominences. The equilibrium of large-scale, 
quasi-static helmet streamers seems to be closelv associated 
with the occurrence of many CMEs te.g.. Uling & Hundhau- 
sen 1986: Hundhausen 1993). and the “disparition 
brusque" phase of prominences normally takes piace in 
changing heimet streamers te.g.. Tandberg-Hanssen 1995). 
Different aspects of these interactions have been treated in 
the literature, but the complexity of the problem generally 
requires serious simplifications in the applied physical 
models. 

It is difficult to construct helmet-streamer solutions ana- 
lytically from MHD equations because both static, mag- 
netically dominated closed regions and flow-dominated, 
open regions occur. Pneuman &. Kopp (1971) first obtained 
a numerical solution for a helmet streamer using an iter- 
ative procedure. Since then, a more general relaxation 
method has been used to obtain quasi-static helmet 
streamers (Steinolfson. Suess. & Wu 1982: Wang et al. 
1993). It has been demonstrated that the helmet streamer is 
a better initial state for simulating looplike coronal mass 
ejections (Steinolfson & Hundhausen 1988; Wang et al. 
1995) as compared to the static initial state because the 
former includes consistently coronai expansion and resem- 
bles the observed features of the preevent corona. 

With a helmet streamer as the initial state, manv attempts 
have been made to simulate the initiation of CMEs, either 
by photospheric shearing i Linker t&; Mikic’. 1995i. magnetic 
flux emergence (Guo ei al. 1991; Steinolfson 1992). or 
thermal energy release in a hcatea corona (Steinolfson & 


Hundhausen 1988), and considerable progress has been 
made. However, little attention has been paid to the internal 
structure of helmet streamers in these numerical simula- 
tions. although there is observational evidence that a helmet 
streamer indeed exhibits complicated magnetic structures, 
e.g., the cavity situated below the streamer dome and the 
quiescent prominence inside the cavity (Pneuman & Orrall 
1986: Low 1994). The magnetic structure in the closed 
region of the streamer is very important because it deter- 
mines both the equilibrium of the streamer and the mag- 
netic free energy in the preevent corona needed to fuel 
CMEs. Aly (1984, 1991) showed that if a force-free magnetic 
field were to be anchored to the surface of the Sun. it cannot 
have an energy in excess of that in the corresponding fully 
open configuration. Low & Smith (1993) interpreted Aly’s 
constraint to mean that, for anchored fields, the field- 
aligned currents cannot store enough energy to fuel the 
CME, and they proposed that magnetic energy in the form 
of detached magnetic fields with cross-field currents may be 
the source of the total mass-ejection energy. They con- 
structed analytically a helmet-streamer solution including a 
magnetic bubble in its closed field region, but without con- 
sidering the plasma flow. Low (1994) further suggested that 
the cavity below the helmet dome may indeed include a 
detached magnetic flux rope associated with the promi- 
nence. In numerical simulations, it is now possible to con- 
struct seif-consistentiy a helmet streamer containing a 
detached magnetic structure via current emepgence from the 
streamer base. 

Recently. Wu. Guo. &. Wang (1995) (hereafter WGW95) 
investigated the dynamics of a helmet streamer by moving 
quasi-statically into the closed field region of the streamer 
field. They showed, with ideal MHD. that if the emerging 
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magnetic bubble is sufficiently small, the resulting configu- 
ration will be a quasi-equilibrium helmet streamer contain- 
ing a detached magnetic bubble in the closed field region. If 
the emerging bubble is large, the solution describes a non- 
equilibrium state, and after a period of slow evolution, the 
streamer-bubble system becomes unstable and breaks up 
dynamically. Their results indicate that a helmet streamer 
containing a detached magnetic structure in its closed field 
region has sufficient magnetic free energy to easily go 
unstable. The quasi-static solution obtained may serve as a 
preevent corona for CME simulations. 

In this paper we shall investigate the dynamic response of 
a helmet streamer to the emergence of a current from below, 
but this time with its magnetic field in a direction opposite 
to the overlying streamer field and include the effect of resis- 
tivity. We do not attempt to model the emergence of a 
specific current, but rather, we are interested in the topo- 
logical change of the magnetic field of the streamer, as a 
means to construct more realistic helmet-streamer solu- 
tions. 


2. NUMERICAL MODEL 
2.1. Governing Equations 

We model the coronal plasma by using the single-fluid, 
resistive MHD equations in spherical coordinates assuming 
axisymmetry. The equations are identical to the equations 
used in WGW95, except for the magnetic induction equa- 
tion and the energy equation which are modified to include 
the effects of resistivity, viz., 

^ = Vx(FxR)-Vx 7 /, (1) 

and 

dT * * * 

* - V(TV) + (2 - y)7W + y -=~ 7 J 2 , (2) 

where the current J = [l/g 0 )V x B, rj is the magnetic resis- 
tivity, and the other symbols have their usual meanings. 
Like Wu. Bao. & Wang (1986). we choose for 7 the follow- 
ing form: 


'o(|J|-|y f |) 2 , \j\>\j t \ 

.0 \J\<\J'\' 


( 3 ) 


where a is a constant which scales with Spitzer’s classical 
resistivity a 0 . and J r is the magnitude of the critical current 
above which the resistivity becomes effective. Physically, 
this corresponds to some microscale instability that could 
enhance the resistivity by several orders of magnitude over 
the classical value when the current intensity exceeds the 
critical value (Song, Wu. & Drver 1989). In our computa- 
tion, we choose a = 6.94 x 10 9 a o and J c = 1.14 x 10~ 6 A 
m . Some comments about the magnetic resistivity in our 
model are presented in § 4. 

The computational domain extends from 1R S to 7R,, 
where R, is the solar radius, and from the pole (6 = 0°) to 
equator {9 = 90°). We use a nonuniform grid in the r- 
direction, Ar, = 0.025r,_ , and a uniform grid in the merid- 
ional direction. Ad = to obtain a grid of 80 x 63 
elements. The MHD equations are discretized by the com- 
bined difference scheme (WGW95), and the method of pro- 
jected normal characteristics (Wu & Wang 1987) are used 
to treat the boundary conditions. At the equator and the 


pole symmetric boundary conditions are used; at the inner 
and outer boundaries characteristic boundary condition 
and linear extrapolation are used, respectively. 

In order to investigate the dynamic response of the 
helmet streamer to the emergence of a current, we first con- 
struct a quasi-equilibrium helmet streamer by using the 
relaxation method (Steinolfson et al. 1982). The parameters 
of the streamer are the same as in WGW95. The plasma 
parameters at the inner boundary are n 0 = 3.2 x 10 8 cm -3 , 
7o = 1.8 x 10 6 K, B 0 = 2.0 G (i.e., the initial value of 
plasma-/? is 1.0 at the equator and 0.25 at the pole). The 
magnetic field lines and velocities of the initial streamer arc 
shown in Figure 1. 

In the second step, we construct a current-carrying mag- 
netic bubble and move it quasi-statically into the closed 
region of the streamer to simulate the current-emergence 
process. To represent the two-dimensional current-carrying 
magnetic bubble, we assume an equilibrium axisymmetrical 
cylinder of radius a containing a volumetric current. Using 
local cylindrical coordinates (r\ 6’, c), and solving the equi- 
librium equations, we obtain the pressure and the magnetic 
field (see WGW95) in the form 

P(r') = Po + Hojlirsa 4 - W r ' 2 + iW ' 3 - -j^r' 4 ) , (4) 

and 

B{r') = - no jo(W ~ ^' 2 )e r , (5) 

where 


■2 18(m - l)p 0 


and p 0 is the plasma pressure at the inner boundary. 

To implement the current-emergence process, wc trans- 
form the expressions for p and B from local cylindrical coor- 
dinates to global spherical coordinates (see WGW95). This 
transformation introduces a small error [V=x=8/(B 0 /R s ) < 
0.005] to. the solenoidal condition of the magnetic field. To 
assure that" the solenoidal condition is satisfied, we have 
employed the iterative divergence-cleaning procedure 
(Ramshaw 1983) in our computation. 

The process of the current emergence is accomplished as 
follows. First, we place the current-carrying magnetic 
bubble below the inner boundary (r = 1R,) with its center at 



Fig. 1. — Magnetic field lines and velocity field of the initial helmet 
streamer. 
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a 7c. ** ~ ?' We thcn movc ix vei T slowl y- with respect to the 
Alfven veloaty, V if into the computational domain at con- 
staat speed (* <$ *). We choose approximately 4 hr for the 
bubble to move completely into the closed region of the 
streamer by adjusting * according to the size of the bubble 
in order to maintain the system in a quasi-equilibrium state 
As soon as the current-carrying structure is inside the closed 
region of the streamer, we let the full set of MHD equations 
evolve and derive the numerical results. It should be noted 
that the lower boundary conditions are also computed at 
each time step m accordance with the method of character- 
istics (Wu & Wang 1987). The numerical results are given in 
the next section. 6 


4.5 R„ while reconnection is still occurring at the flanks of 
the bubble. 

The corresponding relative density contours show that 
the largest density enhancement occurs at the side of the 
loop where the current is large, while the top of the loop has 
a relatively small density enhancement and current. 

To illustrate the later phase of the eruption we have 
plotted field lines and vector velocities at 10 and 20 hr after 
ihe current emergence shown in Figure 4. We see that the 
streamer is disrupted by the passage of the bubble, and that, 
after the disturbance has propagated out of the domain, the 
streamer returns to its initial state atr = 20 hr. 


3. NUMERICAL RESULTS 

In this study we have computed three cases with differen 
sizes of the bubble radius leading to different levels of th< 
SSSy ( ° r ^ urrent ^tensity) contained in these magneti. 
bubbles. Subsequently, we obtain three distinct solutions a< 
shown in Figure 2. representing three different states, viz 
(1) an eruptive, (2) a nonequilibrium. and (3) a quasi- 
equilibrium state. The detailed physical results for the three 
cases are discussed below. 




3.1. Eruptive Case 

. i n u, thiS u Ca f e We choosc = ai 5 *i- the size of th 
bubble which possesses a magnetic energy of 1 57 x 10 3 

ergs, and an initial veloaty * = 14.5 km s" 1 (i.e.. 0.06* 
The evolution of the magnetic field and the velocity as wel 
as of the current and density for t = 4, 5, 6, 7 hr are showi 
in Figure 3. We see that the emerging current quickly dis 
rupts the initial helmet streamer. Because the field lines o 
the emerging field are in a direction opposite to the fieh 
lines of the streamer, a current sheet forms between thi 
upper haff of the bubble and the streamer field. Figures 3a 
3b, and 3c show an amplified view of the lower part of th< 
computational region. It is clear that two x-type neutra 
points form on the flanks of the bubble where the curren' 
(tensity is maximum. Magnetic resistivity causes the fielc 
lines to reconnect, and the velocity at the two ends of the 
-reconnected region shows that the reconnected field i< 
moving away from the reconnection region. The top of the 
erupting current loop reaches a speed of 170 km s~ l ai 



2 ' Evolutlon of thc position of the neutral point at equator 
between preexisting held and the currem-induccd Held in three cases. 


3.2. Nonequilibrium Case 

In this case the bubble size is taken to be 0.07R the 
magnetic energy is 3.41 x 10 29 ergs, and the velocity* is 
6.77 km s 1 (i.e., 0.028*). b 

Figure 5 shows an amplified view of the topology change 
of the magnetic field in the closed region of the streamer for 
this case. We notice that the reconnection process is much 
slower than in case 1. At t = 2 and 4 hr there is no obvious 
reconnection occurring at the bubble sides. However, as 
time elapses the effect of the reconnection becomes clear, in 
that the two x-type neutral points at the flanks of the 
emerging bubble move closer and the emerging current 
slowly dissipates. Finally, the two x-type neutral points 
merge into one and the emerging current degenerates into a 
point. At the same time the current existing above the upper 
half of the emerging bubble expands, forming a new bubble 
structure above the degenerated current. The field direction 
of the new bubble is the same as that of the overlying 
streamer field. The new configuration resembles magnetic 
field conditions in an inverse-polarity prominence (Leroy 
Bommier. & Sahal-Brechot 1984). 

Figure 6 shows the evolution of the magnetic field and 
velocity from 1 to 1R S . Comparing Figures 5 and 6 we 
realize that as the magnetic topology in the closed ration of 
the streamer evolves, the overall streamer structure is 
almost unchanged. After the emerged current degenerates, 
the new configuration will still experience a stage of slow 
evolution. Then, at about r = 24 hr, the whole structure 
becomes unstable and breaks up. forming a very elongated, 
narrow transient. Its relative density contours at 28 hr are 
shown in Figure 7, defining a transient similar to the 
coronal jets described by Burkepile & St. Cyr (1993). After 
the transient has propagated out of the computational 
domain, the streamer structure recovers after roughly 60 hr 
(Fig. 6/). 


3.3. Quasi-equilibrium Case 

In this case we choose a bubble size 0.06 R s giving an 
amount of magnetic energy of 2.51 x 10 29 ergs, and the 
velocity * is 5.80 km s 1 (i.e., 0.024*). From Figure 8 we see 
that the first stage of the evolution is similar to case 2, i.e.. 
the emerging current degenerates into an x-type neutral 
point, and above it the original current sheet expands to 
form a new bubble. However, because of the smaller input 
energy contained in the original emerging current, after a 
stage of slow evolution and restructuring, theifinal configu- 
ration remains quasi-staticaliv in the corona. At this stage, 
magnetic reconnection is still occurring near the neutral 
point below the bubble due to the finite magnetic resistivity. 
Hence the bubble structure slowly dissipates, and the dissi- 
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FIELD lines AND VELOCITY FIELD (T * 10 HOURS) 
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Fig. 4. Evolution oi magnetic field lines and velocity 
T ~ 1 10 hr (top) and T = 20 br {bottom) for 1. 
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4. DISCUSSION 

4.1. Operating Forces 

Tod j scuss the d ? namic processes that cause the bubble- 
destabilization of the streamer, we have plotted in Figure 9 
the three relevant forces, i.e.. pressure gradient, magnetic 
lorce, and gravity, at three locations, as a function of time 
tor case 2. At the lower half of the bubble. Figure 9a , the 
positive pressure-gradient and the magnetic force balance 
gravity and support the bubble. At the upper half of the 
bubble. Figure 9b. and 9c. the pressure-gradient force tends 
to overcome the combined negative magnetic force and 
gravity and stretches the closed field lines of the streamer 
and the force unbalance causes the bubble to move upward. 
On the other hand, in case 3. the three forces achieve 
dynamic equilibrium, and the bubble-streamer system is 
quasi-static. Hence, as pointed out in WGW95. the equi- 
librium or nonequilibnum state of the bubble is determined 
bv the dynamic interaction of the three forces. However, 
unlike the simulation in WGW95. where the whole bubble 


structure emerged from the lower boundary, here the new 
bubble is formed by the reconnection between the emerging 
field and the original streamer field. The magnetic confining 
force, produced by the closed streamer field above the 
bubble is smaller in the present case because part of the 
closed streamer field has reconnected and is no longer tied 
to the solar surface. That is why, in the present case, the 
bubble-streamer system can lose equilibrium more easily 
than in the cases studied in WGW95. We find that the 
equilibrium of the final configuration is very sensitive to the 
size of the original input current. 

One of the noticeable common characteristics in this 
paper and in WGW95 is that the emerged structure, in a 
certain parameter range, will first experience a period of 
slow evolution, then the streamer loaded with this structure 
loses its equilibrium. This study and WGW95 show that 
coronal streamers have the ability to be loaded with extra 
current and mass. The stability of the streamer and loaded 
structure depends on the dynamic interaction of the three 
forces as pointed out above. This slow evolutionary stage 
reflects the restructuring process in the streamer in response 
to the emerged structure and may explain the slow expan- 
sion of streamers prior to their eruption filling & Hundhau- 
sen 1986). 


4.2. Magnetic Topology 

Using the results of the bubble-streamer interactions, we 
have studied the changes of magnetic configurations in the 
closed field region of a streamer. Observations show that 
quiescent prominences are formed under coronal helmet 
streamers fTandberg-Hanssen 1995). and many of these 
prominences fall in the inverse-polarity category (Leroy et 
al. 1984). In cases 2 and 3, we have reproduced the basic 
magnetic configuration of the inverse-polarity prominence 
in the closed region of a streamer (Figures 5, 6, and 8), and 
our primary results indicate that such a configuration can 
easily go unstable. The stability and parameters of this con- 
figuration need further investigation, which may elucidate 
the disparition- brusque phenomenon and its basic relation- 
ship to helmet-streamer insta biiities. 

We have studied the magnetic field development, but a 
word is necessary concerning the electric current 
responsible for the bubble configuration. In a fully three- 
dimensional description the emerging current is connected 
to the photosphere, under which it forms a closed circuit 
through, e.g., the convection zone. In our two-dimensional 
planar configuration we approximate this circuit by con- 
sidering an emerging current forming a closed loop (an azi- 
muthal current) around the solar equator. This is a 
reasonable approximation because the axial dimension (the 
length of a prominence) is much greater than the transverse 
dimension (the width of a prominence), see, e.g., Tandberg- 
Hanssen (1995). We also notice that Low & Hundhausen 
(1995) have pointed out the importance of an azimuthal 
current in the closed field region of streamers, because such 
a current running above a prominence is indispensable for 
the support of the prominence against gravity. 

4.3, Magnetic Resistivity 

The magnetic resistivity in the corona) plasma is very 
small, while fast reconnection processes are needed to 
explain the fast magnetic annihilation process occurring in 
the corona. Numerical simulations have been performed to 
investigate the details of the reconnection in relatively 
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Fig. 5.— Magnetic topology change m the closed field region of the streamer m case 2 


simple configurations te.g.. Forbes &. Priest I9S7). However, 
in our model, the geometry is much more complicated and 
the magnetic reconnection occurs only in a small portion of 
the curved sheet. Hence, m our globai simulation, we must 
assume a reconnection model, (e.g.. the one we have chosen 
in § 2), to simulate the microscaie reconnection process 
occurring in a region that is a few orders of magnitude 
smaller than the mesh size. 

Due to the choice of magnetic resistivity in our model, the 
magnetic Reynolds numoer is about 1000 for case 1 and 400 


for cases 2 and 3. These Reynolds numbers are much 
smaller than the magnetic Reynolds number in the real 
undisturbed corona. The problem of unavoidable increase 
of resistivity in numerical simulation of magnetic reconnec- 
tion in solar plasmas is widely recognized and till now, no 
simulation can accommodate the actual magnetic Reynolds 
number under coronal conditions. We have performed the 
simulations using different magnetic resistivity values. We 
tound that different values of tj only affect the solution 
quantitatively, while the basic scenario for the streamer 
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